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Abstract. We present a summary report on a detailed study of the mass 
and dynamical state of the Corona Borealis supercluster using the largest 
redshift survey of this object to date. Extensive A-body experiments 
are performed on simulated superclusters in both critically-flat and low- 
density cosmological models to test the robustness and errors of the mass 
estimators used on the data. From the redshifts and dynamics of 528 
galaxies in the supercluster, the mass of Corona Borealis is estimated to 
be at least 3 x 10^® h~^ Mq, and the mass-to-light ratio is 560 h {M/L)q, 
yielding a matter density parameter of Dm ~ 0.4. We also discuss the 
prospects for mapping the mass distribution on supercluster scales with 
gravitational leasing. 


1. Introduction 

Among the most remarkable scientific discoveries of the past decades is the 
realization that most of the mass in the Universe may reside in dark matter. 
Although ongoing direct searches for dark matter in particle detectors have not 
yet revealed positive results, compelling astrophysical evidence for its existence 
is now abundant. Moreover, various observations now indicate that Nature’s 
ability to hide dark matter increases with distance scales. A useful measure of 
this “hiding power” is M/L, the ratio of mass to light of an object expressed 
in the solar unit {M/L)q. A larger mass-to-light ratio therefore indicates the 
presence of more dark matter for the observed amount of light. Current mea¬ 
surements give a mass-to-light ratio of about 5 in the solar neighborhood, 10 in 
the cores of elliptical galaxies, 100 for the Local group, and as high as 200-300 
for clusters of galaxies (e.g., Carlberg et al. 1996). Nature’s appetite appears 
to be more voracious on grander scales, gobbling up large quantities of matter 
while exhibiting little sign of digestion. 
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The motivation for our current study is to obtain measurements of galaxy 
and matter distributions on the little-explored supercluster scales, and to inves¬ 
tigate whether the mass-to-light ratio continues to rise beyond the much-studied 
cluster scales of a few Mpc. Superclusters are among the largest structures ever 
registered on maps of the nearby galaxies. Each of these immense clusters of 
clusters consists of possibly 10,000 or more galaxies in a region of space tens 
of Mpc across. Our study is focused on the most prominent example of su¬ 
perclustering in the northern sky - the Corona Borealis supercluster, centered 
at right ascension 15^20™ and declination -|-30°. This extraordinary cloud of 
galaxies was already noted by Shane and Wirtanen in 1954 and Abell in 1958. 
The Corona Borealis supercluster is now known to consist of seven rich Abell 
clusters at z ~ 0.07 and numerous galaxies in the intra-cluster regions. Its core 
covers a 6° x 6° region of the sky, corresponding to nearly 20 x 20 Mpc in physical 
size. Whether the supercluster extends beyond this scale is not yet clear; deeper 
and wider surveys of the sky will be required to answer this question. 

The previous redshift survey of the Corona Borealis supercluster was fo¬ 
cused on the densest regions surrounding the cores of six individual Abell clusters 
A2061, A2065, A2067, A2079, A2089, and A2092 (Postman, Geller, & Huchra 
1988). From the dynamics of more than 150 cluster galaxies, the sum of the 
masses within the central lh~^ Mpc of each of the six clusters was estimated to 
be 2.4 X 10^® AIq. The mean mass-to-light ratio in the R band for the same 
region was about 250. If the same M/L is assumed to extend to regions beyond 
1 Mpc, the supercluster mass is then 8.2 x IO^^/i^^Mq. A higher M/L ratio 
would further increase the mass. 

A more extensive redshift survey of Corona Borealis has since been com¬ 
pleted with the 176-fiber Norris Spectrograph at the Palomar 200-inch telescope 
(Small et al. 1997ab; 1998). This survey substantially increases the number of 
galaxy redshifts both in the Abell clusters and in the intra-cluster regions. It 
allows for a more accurate determination of the structure of the supercluster and 
a more reliable estimate of its mass. A total of 32 fields, each with a 20-arcmin 
diameter field-of-view, are successfully observed, yielding redshifts for 1491 ex- 
tragalactic objects. With 163 redshifts from the literature, the entire survey 
contains 1654 redshifts. A total of 528 galaxies actually lie in the redshift range 
of the supercluster, 0.06 < z < 0.09, and forms the sample for the study reported 
here on the supercluster’s mass and structure. Figure 1 shows the redshift and 
right-ascension pie diagrams in three declination slices for all galaxies in our 
survey with z < 0.15. The supercluster is well delimited along the line of sight 
by foreground and background voids. A background supercluster at z ^ 0.11 
including Abell 2069 is also visible. The foreground galaxies at z ~ 0.03 are 
part of the “Great Wall” identified by Geller and Huchra (1989). 


2. Testing Mass Estimators 

The main difficulty in estimating the mass of a supercluster is that unlike a star 
or galaxy cluster, a supercluster is so large that the galaxies and clusters residing 
in it have not had sufficient time to interact and randomize their motion. It is 
therefore unclear that traditional estimators such as the virial theorem are valid 
for a supercluster. In order to address this problem, we have chosen to test the 
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Figure 1. Redshift and right-ascension distribution of all galaxies 
with redshift z < 0.15 in our Corona Borealis survey. The declination 
range covered by the survey (27.5° < 5 < 32.5°) is shown in three 
slices. The thick lines mark the redshift extent of the Corona Borealis 
supercluster, while the dotted lines mark the background supercluster. 
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performance of the virial and other mass estimators on computer-simulated su¬ 
perclusters drawn from large A/^-body simulations of cosmic structure formation. 
The simulated superclusters, which are in general quite spatially anisotropic, also 
enable us to assess the effects of the non-uniform sampling in our observations 
on our mass estimates. 

To anticipate the possibility that supercluster dynamics has a systematic 
dependence on the cosmological density parameter, we have simulated both 
critically-flat and low-density models. The = 1 rnodel is the standard cold 
dark matter (CDM) with a Hubble constant of h = 0.5 and a normalization of 
(Tg = 0.7 for the rms mass fluctuations in spheres of radius 8h~^ Mpc. The 
low-density model is CDM with = 0.3, a cosmological constant = 0.7, 
and h = 0.75, normalized to the 4-year COBE quadrupole Qrms-PS = 18 
(Gorski et al. 1996). The corresponding ag is 1.2. We performed both large- 
box simulations (640 Mpc a side) with random Gaussian initial conditions and 
small-box simulations (160 Mpc a side) that are constrained to produce objects 
with an overdensity of roughly 5. A comoving Plummer force softening length 
of 160 kpc is used in the particle-particle particle-mesh (P^M) force calculation 
for all simulations. 

Figure 2 shows the dark matter distribution in one simulated supercluster. 
Each side of the box shown is 30 Mpc, and the total mass is 1.1 x 10^® Mq- 
This object is extracted from a large cosmological simulation of the CDM model 
that traced the motion of 17 million particles in a (640 Mpc)^ box. The average 
overdensity in this region is 5p/p = 4.7. We identified 15 other supercluster¬ 
like objects that emerged in our cosmological simulations of both = 1 and 

= 0.3 models. To minimize a potential selection bias, only two criteria 
are applied in identifying these objects. First, the candidates must consist of 
multiple dark matter halos that have not yet merged into one dominant smooth 
halo. This feature is strongly suggested by the observed complex substructure 
(e.g., the Abell clusters) of the Corona Borealis. The only other criterion is 
that the objects have an average overdensity of ~ 5 in a volume of ~ 30^ 
Mpc^, corresponding to the least evolved dynamical state of Corona Borealis 
(see below). This is a conservative choice since systems with higher overdensities 
are closer to virialization, and the mass estimators should only work better. 

For the virial estimator, the mass is calculated from 


where a is the line-of-sight velocity dispersion and (1 /vp) ^ is the mean harmonic 
projected separation. A simple estimator for (l/rp)“^ is given by 
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where 6ij is the angular separation of galaxies i and j, D is the radial distance, 
and N is the total number of galaxies. This estimator, however, is very sensitive 
to close pairs and is thus quite noisy, especially for systems which have not been 
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Figure 2. An example of the simulated superclusters used to test the 
mass estimators applied on the actual data. This complex structure 
consists of at least 7 cluster-type subclumps, each with a mass exceed¬ 
ing 2 X 10^^ Mq. Each side of the cubic box is 30 Mpc, and the total 
mass inside the box is 1.1 x 10^® Mq. 


uniformly sampled spatially. We adopt an alternative estimator introduced by 
Carlberg et al. (1996) that is less sensitive to irregular sampling and close pairs: 
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where rj and rj are the projected radii of objects i and j, K{k) is the complete 
elliptic integral of the first kind in Legendre’s notation, and kfj = 4rjrj/(rj+rj)^. 
Although this “ringwise” estimator was originally developed for systems such as 
galaxy clusters with near circular symmetry on the sky, we find this estimator to 
give less biased values of (l/rp)“^ than the straightforward sum in equation (^). 

We have also tested the projected mass estimator (Bahcall &: Tremaine 
1981) given by 

= ( 4 ) 

i 

It is designed to give equal weights to particles at all distances (if oc 1/r on 
the average), but the estimate depends on the mean eccentricity of the particle 
orbits parameterized by /pm- It can be shown that /pm = 32/n for isotropic 
orbits and 64/7r for radial orbits, independent of the mass distribution (Heisler, 
Tremaine, & Bahcall 1985). We have chosen to use /pm = 32 jn since this yields 
the smallest masses. 

Table 1 summarizes our test statistics on the 16 simulated superclusters 
in both Qrn = 1 and 0.3 models. Each supercluster is projected onto a two- 
dimensional sky, and the masses are computed from eqs. (1) and (4) using all 
the particles in the supercluster. Overall, we find that the virial theorem works 
remarkably well at recovering the true mass Mt of the simulated superclusters. 
An important result is that the estimators do not perform worse in the 0^ = 0.3 
model as one may have expected for a low-density model. 

To test the effects of irregular sampling in our redshift survey, we have 
also conducted simulated observations of these model superclusters with strate¬ 
gies similar to those used during the actual observations. Each supercluster 
is projected onto a two-dimensional sky and portions of it are viewed through 
randomly-placed fictitious observing fields. Instead of using all particles, a frac¬ 
tion of them in each field are randomly rejected so that the total number of 
particles used in the simulated observations was roughly 500, comparable to the 
number of galaxies with measured redshifts in Corona Borealis. We find the 
non-uniform sampling to systematically underestimate (l/rp)“^ but does not 
significantly affect the velocity dispersion. Overall, the virial estimator underes¬ 
timates the true mass by 31% in the Clm = 1 rnodel and by 5% in the = 0.3 
model. Based on this result and Table 1, we conclude that the virial and pro¬ 
jected mass estimators can be reliably applied to superclusters even though they 
are not relaxed systems. With our irregular sampling, the estimators may un¬ 
derestimate the true mass by up to 30%. 

An interesting question to ask is why the virial (or similar) mass estimator 
works so well for unrelaxed objects with prominent substructure such as in 
Figure 2. An important clue is that our 16 simulated superclusters all turned 
out to be gravitationally bound even though this was not a selection criterion. 
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Table 1. Test results of mass estimators 


Model < My/Mx < Mp/Mj' 

= 1 0.94 ± 0.24 0.76 ± 0.24 

Qm = 0.3 0.99 ± 0.21 1.04 ± 0.52 


“Ratio of virial mass Mv of eq. (1) to true mass. 
*’Ratio of projected mass Mp of eq. (4) to true mass. 


How well does the virial estimator work for a bound system? For a random 
collection of particles, the virial theorem clearly can yield a mass that is wrong 
by an arbitrary factor. For a bound system with ii' — | VF| < 0 ( where K and W 
are the kinetic and potential energy), however, it is easy to show that the virial 
mass My can never exceed twice the system’s true mass Mp- (Recall a virialized 
system has K = \W\/2.) In fact. My approaches 2Mt only for a marginally 
bound system with K ^ \W\. For systems with |VF|/2 < K < \W\, one gets 
Mp < My < 2Mp; while for systems with K < |VF|/2, the virial estimator 
underestimates the mass: My < Mp. The virial theorem can therefore recover 
the true mass of a bound system to within a factor of 2 in most cases. We indeed 
find this statement to be true for all 16 simulated superclusters. 

Although our Wbody experiments strongly suggest that superclusters are 
generally bound objects, we provide an additional argument here supporting the 
possibility that the Corona Borealis supercluster is gravitationally bound. The 
spherical tophat model for gravitational collapse indicates that an outer shell 
of a density perturbation is bound if the mean overdensity within the shell is 
greater than — 1)/(1 + z) (e.g., Peebles 1980). Taking 0.3 to be a rough 

lower limit on a supercluster at z = 0.07 must be overdense by more than 

a factor of 2.1 for it to be bound. In comparison, the observed overdensity 
in galaxy counts in Corona Borealis relative to the field is 6 ~ 7/, where the 
dimensionless parameter /, defined by / = c5z/H^dr, measures the relative 
elongation in redshift space 5z and real space 5r along the line of sight. The 
minimal value of / is unity, corresponding to zero peculiar velocities, whereas a 
larger / indicates a stronger “finger-of-god” effect due to galaxy peculiar motion. 
Although we cannot completely rule out the / = 1 case, it is very difficult to 
arrange an object that is as immense and prominent as the Corona Borealis while 
at the same time is freely expanding with the Hubble flow without exhibiting 
any peculiar velocities. The overdensity in galaxy counts in the Corona Borealis 
is therefore at least 7, but more likely to be significantly above this value. The 
Corona Borealis supercluster is therefore bound unless galaxies are significantly 
biased relative to mass by more than a factor of 3.5. 

3. Mass and M/L of Corona Borealis 

Given the robustness of the mass estimators found in our numerical tests de¬ 
scribed in the previous section, we now proceed to apply eqs. (l)-(4) to the actual 
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data from the Corona Borealis survey. Using the techniques of Beers, Flynn, 
& Gebhardt (1990), we estimate the centroid velocity of the Corona Borealis 
supercluster to be cz = 22A20'^\^l km s ^ and the dispersion to be a = 1929~^q^ 
km s“^ in the cluster rest frame. Using the ringwise estimator in eq. (3), we find 
(l/rp)“^ = 4.6 h~^ Mpc. The virial mass estimator (eq. (1)) then gives a mass 
of 3.8 X 10^® h~^ Mq for the Corona Borealis supercluster. The projected mass 
estimator in eq. (4) yields a similar value, 4.2 x 10^® h~^ Mq. Based on the tests 
described in the previous section, we conclude that a secure lower bound to the 
mass of the Corona Borealis supercluster is 3 x 10^® h~^ Mq. 

The supercluster luminosity function for M{Bab) < —16.3 + 5 log h mag in 
the AB-normalized B band is presented in Small et al. (1997b). Integration of 
the luminosity function yields a mean luminosity density of pl{Bab) = 1.9 x 
10® h Lq Mpc“^ in the supercluster. Taking the solid angle of the survey to 
be 0.0076 sr (= 25 deg^) and the redshift limits of the supercluster to be at 
z = 0.06 and 0.09, the volume of the region surveyed is 2.8 x 10*^ h~^ Mpc^. 
The M/L ratio of the supercluster in the Bab band is thus 560 /i(M/L)q. 
For a local BAB-band luminosity density of (1.8 ± 0.2) x 10®/iLq Mpc“^ for 
M{Bab) < —16 + 51og/i mag (Small et al. 1998 and references therein), the 
corresponding critical M/L ratio to close the universe is 1550 ± 170 /i(M/L)q. 
We therefore obtain Llm = 0.36 for the density parameter on supercluster scales 
of ~ 20 h~^ Mpc. 


4. Future Prospects 

Our determination of M/L and on supercluster scales can be strengthened by 
measurements of other superclusters with techniques similar to those described 
above, or by independent measurements with alternative mass-estimating meth¬ 
ods. The former will be achieved when the ongoing large-area redshift surveys 
(such as the 2dF and Sloan surveys) generate sufficient data. For the latter, a 
very powerful technique is gravitational leasing, which probes the mass profile 
of a foreground object by the gravitational distortion it induces on the images 
of background objects. A galaxy cluster, for example, induces a correlated el¬ 
liptical distortion in the images of background galaxies. Measurements of such 
statistical signals over many galaxies have been used successfully to map the 
surface mass densities of individual clusters (e.g., Tyson et al. 1990; Kaiser & 
Squires 1993; Mellier et al. 1996 and references therein). Extending the gravi¬ 
tational lensing technique to superclusters will be a promising alternative probe 
of the dark mass distribution on such large length scales. 

Figure 3 is an illustration of the lensing signal expected in a superclustering 
region. The left panel shows the particle distribution in a different projection 
of the same simulated supercluster in Figure 2. Each bar in the right panel 
shows the local strength and direction of the theoretical elliptical distortion in 
the shapes of the background galaxies due to lensing by the supercluster. More 
explicitly, if the two-dimensional gravitational potential is denoted by , and 
7 i = (5^ — dy)ip/2 and 72 = dxdyip, then the magnitude of each bar is propor¬ 
tional to the shear 7 = ( 7 ^ -|- 7 ^)^'^^ that measures the anisotropic stretching 
of the images. The angle (5 of each bar relative to the horizontal axis obeys 




Figure 3. Left; A different projection of the simulated supercluster 
shown in Figure 2. Each side of the box is 30 Mpc, and the total 
mass in the box is 1.1 x 10^® Mq. Right; The local magnitude and 
orientation of the shear field induced by the supercluster. 


tan(2/3) = 72 / 71 . Comparison of the two figures shows that the amplitude of 
the distortion on the right is clearly correlated with the projected mass density 
on the left, and the preferential tangential alignment of the shear orientation is 
evident. A promising feature is that the shear amplitudes do not drop appre¬ 
ciably along the ridges connecting individual cluster-like clumps. Measurements 
of the distortion in the images of background galaxies can therefore in principle 
be used to construct a surface mass density map in a supercluster. This offers 
a powerful alternative way to probe the distribution of dark matter on scales of 
tens of Mpc. 
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